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Abstract 

This paper deals with the convolution powers of the characteristic func- 
tion of [0, 1], X[o.i] a-nd its function-derivatives. The importance that such 
convolution products have can be seen, for an instance, at [2] where there 
is the need to find the best difFerentiable splines as regularization tool to 
be used to expand functions or in the spectral analysis of signals. Another 
simple application of these convolution powers is in the construction of a 
partition of the unity of a very high class of differentiability. Here, some 
properties of these convolution powers have been obtained which easily led 
to write the algorithm to produce convolution powers of X[o,i] a-nd their 
function-derivatives. This algorithm was written in calc, the program 
published under GPL that is free to download. 
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1 Introduction 

The n-th convolution power of the characteristic function of an interval is 
an (n — l)-b-splines, this is the subject of this paper. This first section will be 
a general description of the work done, the notation used, finishing with three 
lemmas and a theorem, which are the main tools used in this work. 

Then we proceed in section [2] with the calculus of the powers 2,3,4 in a way 
to uncover the recursion which will be finally clear in the calculus of the third 
power. The forth power will be a preparation for the final theorem (algorithm) 
which will be presented at the section [3l The section [2] is, thus, preparatory for 
the proof at section [Sj and was thought to permit the reader to construct the 
algorithm for himself, so the reader can skip to the section [3l if he find that this 
subject is familiar and if he want only to see the main result. 

In the section [4] it will be presented a computer program which is a machine 
translation of the theorem stated at section [31 Finally, the last section, contains 
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some applications suggested for the material in this paper and a characterization 
of a general space of univariated splines. 

We have used calc P], a computer language specialized to construct math- 
ematical algorithms with infinite precision this meaning that it can deals with 
arbitrarily long integers which is fundamental to this work as it will be necessary 
to use factorial for any value of n in the calculus of the convolution powers. In 
section |4] you can see a graphic output of the program for n = 30. 

1.1 Some history 

Splines are pieces of polynomials pasted together in conditions to make a 
continuous, and even more, a highly differentiable function. In the last section 
you may find a definition for a class of splines. It is a simple idea which turned 
into a very powerful tool in such diverse areas as geometric design or approxima- 
tion theory. In the background there is a famous theorem of Weierstrass stating 
that any continuous function can be arbitrarily approximated by polynomials. 
There was a huge work to make this simple and easy, and we can consider a bad 
example the Taylor polynomials, which make a very good approximation of a 
function in a sole point but which is good enough to be used to approximate pe- 
riodic functions like sine. The bad thing with classical polynomial approximation 
is that computers do not like it that much because it needs a very high degree 
to get something useful. You can remember Lagrange polynomial interpolation 
which does a good job if the points to be linked together are not jumping to 
much and the critical situation is at the boundary of the interval which defines 
the region where the points have been selected. This is a particular problem 
to classical polynomial approximation because just around the boundary of the 
region containing its roots a polynomial will behave wildly. If you have selected 
n points the corresponding Lagrange polynomial is of degree n — I. 

Then splines came out to be the solution. With some pieces of first degree 
polynomials you may have a very good approximation of data collected from a 
sensor, provided that you are only interested in the total amount of the phe- 
nomenon, its integral. Well this was the trick, paste together several piece of 
polynomials instead of a unique monolithic polynomial as a Lagrange polyno- 
mials. The piecewise first degree polynomials are 1-splincs, known as polygonal 
lines. 

Switching to polynomials of second degree we can already have a better ap- 
proximation and the famous algorithm of Runge-Kutta does exactly this for 
solving differential equations, though not with splines. But degree two polyno- 
mials have a big lack of flexibility since they do not provide a change of convexity 
in each piece, they are 2-splines. So, the next step is to paste together pieces 
of polynomials of third degree and this has proved to be a good solution for a 
large number of problems needing polynomial approximation: 3-splines. 

At this point it should pop up a definition for this little thing called splines 
that started to be used in the seventies and today is a very big network of 
research in Mathematics. To be an n-splines, the pieces of polynomial should 
be of degree less or equal to ??, there can be a straight line infiltrated in the 
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business. The second condition is that the resulting function has to be of class 

It is simple like that but the main literature about splines, for example, the 
historic book of Larry Schumaker [5] , is not easy to be read though it is at the 
beginning of the history. This is the point with this paper, it is starting by the 
end, not from the beginning, and there are some powerful tools that had been 
added to the history to make splines easy. 

1.2 The contribution of this paper 

In chapter 6 of [1], the author has shown that the n-th convolution power 
of the characteristic function of the interval [0, 1] is an (n — l)-splines having 
as support the interval [0,n] and he starts by showing that the convolution 
of two characteristic functions results in a 1-splines. Here you are the basic 
tool to build splines: the convolution. Convolution is a "multiplication", it 
satisfies the distributive law relatively to the sum, hence a sum of characteristic 
functions of intervals can be multiplied by a selected characteristic function of 
an interval, using convolution, to produce a 1-splines. Iterating this process 
you will reach an n-splines for the desired n. But this take a lot of time and 
computer power and the alternative is to construct a splines-bases to generate 
an appropriate space of splines and this is the tool that this paper is going to 
produce: an approximate unity to build up bases for vector spaces of splines. 
This is an well known application that will be only mentioned at the last section 
as complement to the main result of this paper. These convolution powers can 
be easily be transformed into a sharp approximate unity with respect to the 
convolution product and this will be mentioned again at section [S] 

Here we have a slightly different methodology from the one at [1], besides 
presenting four properties of / = (x[o,i])"i which arc the lemmas and the the- 
orem at section [21 the work is directed to a recursive program that is able to 
produce the whole formula not only of / = (x[o,i])" but for all its function- 
derivatives, having as input the power n which is needed. This is the main 
result of the paper, theorem HI in the section [3] followed by the correspondent 
algorithm, a computer program, in sectional The program produces functions 
for gnuplot [7], which may be run to show the graphics of / and is able to write 
the matrix of the coefficients of / and all its function-derivatives. 

There are recent works, see [2] for example, which show the importance of 
convolutions in "spectral analysis" mainly in the context of wavelets. This is an 
example of the importance of having a simple way to present convolution powers 
and a good algorithm to produce them. This paper offers that contribution in 
this matter in the particular case of convolution of characteristic functions of 
intervals of the line. 

The whole construction is done step by step, in a way to permit the reader 
to construct by himself the method. 

Our main motivation with this paper is the construction of a tool to be used 
in the analysis of differential operators and in the spectral analysis of wavelet 
and we think this is underway. 
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2 The convolution powers of x 

In this section, first, we introduce the notation and basic results necessary 
to proof the main theorem of the paper, following this we are going to calculate 
the powers 2, 3, 4, in order to formulate the general result. 



2.1 Basic results 

Here you are the basic notation that will be stated together with three Lem- 
mas which are the underground tools for the theorems of the paper. 

The characteristic function of the interval [0, 1] will be referccd to using one 
of the following notations: 

Xo = X[o,i] = X (1) 

It will be used the index when necessary to indicate a translation of a G R, 
then 

Xa(a;) = Xo(a; - a) = X(a; - a); (2) 

The use power notation for "convolution power" has no risk of possible 
confusion with the "arithmetic power" , since the context will indicate what is 
the meaning in the few occasions in which the "arithmetic power" will be used. 
Hence X = X^ is the first power, as usual, and we can add a sense to as the 
constant function. 

The symbol / = x" '^iU be used to simplify the notation of the n-th con- 
volution power, which is the main goal of this paper, hence / do represent the 
n-th convolution power of Xi consistently, in the whole paper. 

The second convolution power of the characteristic function is the triangle 
function, a 1-splines, whose derivative is the the difference of two the translations 
of the characteristic function: xo ~ Xi ■ This is a consequence of the following 
lemmas. 



Lemma 1 (Convolution and derivative) 



4-{f *9) = 3^*5 = /*^ 
ax ax ax 



Lemma 2 (Distribution derivative of characteristic function) 



Lemma 



— = do - Oi 
ax 

3 (Convolution with 5a and translation) 

{5a*f){x)^f{x-a) 



Using these lemmas it can be proved: 
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Theorem 1 (Derivative of the second convolution power) // 



then 



f = Xo - Xi 



In lemma [31 when the translation parameter is a = 0, we can see that the 
Dirac delta measure is the unity of convolution as an element of an extension 
ring of the convolution ring of functions which does not have a unity. 



2.2 The second convolution power 

These three properties stated as lemmas produce an easy way to get the 
second convolution power of x, the above mentioned triangle function, (compare 
with the heavy work at [1] chapter 6]. Moreover, this method will be used 
all along the paper, repeatedly, the main reason being that we shall have a 
succession of integrals whose domains of integration are disjoints, they coincide, 
in fact, exactly on a single point. 



/' = (x * X) = X * X = X * (<^o - (^i) 
/' = (X * X)' = X * (^o - Xo * Si 
f = (x * x)' = X - Xi 

fix) = X * X{x) = x\x) = / {x{t) - xi(t)) dt 



fix)^xHx) 



fix) = x'^ix) = < 



fix) = x^{x) 








X <0 







X 


X e [0, 1] 


Jdt 






X 


X e [1,2] 


«-/ 
1 


x>2 





' x<0 





X e [0, 1] 


X 


[1,2] 


1 - 


^ x>2 





X <0 





X e [0, 1] 


X 


X e [1,2] 


l-( 


x>2 






(3) 
(4) 
(5) 

(6) 



(7) 



(8) 



(9) 



and now it is easy to see that x^ is the so called triangle function centered at 
1 with support [0,2] and height 1 at 1. You may observe the presence of the 
symbol "a", at equation (7), which is the value of the integral on [0, 1] and is 
the initial condition of the integral on [1,2]. This trend will be repeated further 
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and further for all convolution powers. The symbols "a, 6, c, . . ." , will be used 
temporally, before their values can be expressed. 

The triangle function can be defined as x^(a;) = xx + {'2 — x)xi- This expres- 
sion is not useful in computations, but it will be of good help in the mathematical 
formulation, in the logic preparation which will produce the algorithm. 

Observe that the (distribution) derivative of x is '^o^'^i where Sa is the Dirac 
measure concentrated at the point a. The fact that the "convolution with the 
Dirac measure is a translation of the other convolution factor" has been used in 
equation (5). 

An easy and beautiful theorem, extending theorem [TJ can be stated, imme- 
diately. 



Theorem 2 (Derivative of /) the n-th convolution power 

the derivative of a convolution power f of x is the difference of translates of the 
previous power. The proof is a direct application of Lemmas [TJ [2] and [31 This 
theorem will be used several times in the sequence. 

A notation will be introduced to simplify the expression of /. Put Pi{x) = x 
then the equation (9) reads: 

.T < 

f(x)'Y'(x)-i -^^[O'l] ^l(^) (11) 

x>2 

Observe that the second convolution power has been defined in terms of a 
linear combination of a polynomial of degree one and its translates (the plural 
will be valid in a near future. . . ), this highlights that /, the second convolution 
power, is 1-b-splines as particular case of the first sentence of the paper, that 
"n-th convolution power is an (n — l)-b-splines" . It will be possible to go further 
using this kind of polynomials and the same will be done with the next power 
so to make clear the recursion of the process. 



2.3 The third convolution power 

Before going ahead let us stress that the whole secret is contained in first 
calculating the derivatives of the next power which we already have, thinking 
that, to calculate the next power first we must have the previous one. . . So let's 
calculate the derivatives of the third convolution power to have: 



£,X^{x) = x(x) - 2x(x - 1) + x{x - 2); 



(12) 
(13) 
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where we can see the binomial coefBcients: 1,-2,1 as usual when calculating 
derivatives of products and, here, there is a "difference" in the "game" . 
This suggests, by the way, another beautiful and easy theorem: 



Theorem 3 (The n — 1 derivative) of the n convolution power 



1. The n-th convolution power, f ~ is an (?? — \)-b-splines of compact 
support [0,n], hence f has n — 1 function-derivatives of whose n — 2 are 
continuous. 

2. The derivative of order n — 1 of the n-th convolution power f , which is 
not continuous, is given by the linear combination of translates of x 



j(n-l) f 



The coefficients of this linear combination are the elements of the line of order 
n — 1 of the "alternate Pascal Triangle", these elements are the coefficients of 
{x — 1)"~^, (the line of order zero has only one element, 1). 

The proof comes by recursion with repeated application of the lemmas [U [5J 
[3] and theorem [TJ 

The calculus of the iterated integral of /" will be easy as it is again a 
succession of integrals whose integration domains arc disjoints (coincide on a 
set of measure zero), one point set. 

Hence the third convolution power is twice the integral of its second deriva- 
tive starting at the initial condition a; = 0. We know in advance that the support 
of this convolution is [0,3], from [1], and that there will be different equations 
each time we go over an integer boundary. The initial condition of each integral 
is the total value of the previous integral to make a continuous function. This 
will be represented by the symbols a, 6, temporally. 

In condition to obtain the recursion formula, it is necessary to introduce here 
a rather complicated notation, but intuitive, indeed. It will protect the parts 
of the expression which will pop up. Observe that the second derivative of the 
third convolution power is not continuous as a linear combination of translations 
of X with coefficients being the alternate binomial numbers, as in the powers of 
(x-1). So let the following numbers be defined: 

r Bin{n,k) = {-lf{l) 

I aoo = Bin{2,0); 

\ am = Bm{2,l); 

y ao2 = Bin{2,2); 

With this notation it may be written 



(14) 
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X' = fix) 



X <Q 

X <l 

X <2 
x<3 
x>3 





aoixi^ - 1) 
ao2Xix - 2) 




(15) 



Integrating twice this second derivative we shall obtain the formulation of /. 
At this point, / is the second convolution power of x- A new set of coefficients 
will be defined in the run and, for a while, the convention, state previously, to 
use the symbols a, b as the initial conditions for the second and third integrals, 
will be used. 



fix) = lo (aoox(i) + daixit - 1) + «02x(i - 2)) dt 
' T- ^ n n 

fix) = < 



a; < 

.T < 1 

x <2 

X <3 
x>3 



fix) = 



X <0 

X <l 

X <2 
x<3 
x>3 



fix) = 



x<0 
x <1 
< x<2 
X < 3 
X > 3 





aoox 

a + aoi xit - ^)dt 
b + ao2 % xit - 2)dt 




aoQX 

a + aoi xit)dt 
b + ao2 Jq^^ xit)dt 




aoQX 

ooo + aoiix - 1) 

ooo + floi + 002(2; - 2) 





fix) 



X <0 

X <1 
< x<2 

a; < 3 

a; > 3 





awx 

an + aoi{x - 1) 
ai2 + ao2(x - 2) 




oio — ooo; 

aii= aio-Pi(l); 

ai2 = an + aoiPi(l); 



(16) 
(17) 



(18) 



(19) 



(20) 



(21) 



The equation (21) can be written a bit different in a way to show that 
the sum of certain indexes is constant in each row. This rule is valid with two 
exceptions^ which are meaningless to the definition of the algorithm: 



^though that, at this point, the number of exceptions is very high relatively to the number 
of cases, 1, for which the rule is in force . . . 
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1. The first row of the matrix aoj = Bin(n,j); 

2. Tlie two first elements in in each row 

O-iO = flOO; 111 = fljo = floo 

aio = aooA(l); 
an = 010/^0(1); 
ai2 = aiiPo(l) + aoi-Pi(l); 

In equation (22) it is being used the family of polynomials 

Pkix) = x^/k\- 



(22) 



(23) 



In the equations to define the coefficients aij we have in each row a sum of 
products ay Pfc(l) and the sum i + fc is the order of the row when j > 2 (the 
above mentioned exception). 

To obtain a new row Oy- , we have to use all the possibilities in "i + k =order 
of the row" . This law will be valid in the sequence. 

Calculating the next integral we shall have: 



a; < 

X 

X < 1 oio / tdt 



X < 2 a + fli 



X X 

111 J dt + gqi J{t — l)dt 

1 1 

X X 

x<3 b + ai2jdt + ao2j{t-2)dt 

2 2 

X > 3 



X <0 

X <1 010/2(2;) 

X <2 



x—i x—i 

010^2(1) + On / dt + aoi J tdt 


x-2 x-2 

X < 3 b + ai2 J dt + ao2 / tdt 


.T > 3 

a; < 

a; < 1 020^2(2;) 

a; < 2 020^2(1) + anPi{x - 1) + aaiP2{x - 1) 

a;<3 fe + ai2Fi(a;-2)+ao2P2(a;-2) 

a; > 3 

V 



< X 

X 



X <0 

X <l 020^2(2;) 

< 2 a2o-P2(l) + anPi{x - 1) + 001^2(2; - Ij 

< 3 020^2(1) + aii-Pi(l) + 001^^2(1) + ai2Pi 
> 3 



a; > 3 



(a; - 2) + 002^2(2; - 
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' a; < 

a; < 1 020^2(2;) 

f{x) = < x<2 a2iP2{l) + anPi{x 

X <3 022 + 012^1 (.T - 2) - 

x>3 



- 1) + 001^2(2; -1) 
002^2(2: - 2) 



Some remarks here will make things clear in the next step. In equation 
(24) the symbols a, b have been used as explained before, and in equation (25) 
the limits of integration have been changed accordingly with the interval of 
definition, and the liberty of translation invariance of Lebesgue measure has 

.T-2 

been used to change the limits of J dt in conditions to have the same setting 



in both integrals. This leads to a translation of the polynomial P,;, naturally. 
This will be further commented in a moment. The symbol a has been evaluated 
as /(I) and the symbol b has been evaluated as 

b = /(2) = 021^2(1) + aiiPi(l) + 001^2(1) 

the value of the primitive over [1,2] evaluated at point x = 2. On the interval 
[2, 3] it can be written 022^0(2; — 2) to enhance the property 

amjPiix — k)]m + i = j 

where j is the order of the row. The new set of coefficients will be defined 
explicitly: 

{020 = flio = floo; 
021 = 020^2(1); (30) 
022 =a2iPo(l)+aiiA(l)+ 001^2(1); 

and remember the exception to the rule that the first two coefficients in each 
row have a particular law. With these definitions / is 



a; < 

X <1 

X <2 
x<3 
x>3 





020-^2(2;); 
a2iPaix - 1) 
a22Poix - 2) 




anPi(x-l)+aoiP2(a;-l) (31) 
ai2Piix - 2) + 002^2(0; - 2) 



a linear combination of the elements of last row of the matrix a[i, j] with trans- 
lations of all the polynomials P^; fc = 0, . . . , n, where n + 1 is the desired convo- 
lution power. The previous raws of this matrix give the derivatives of / but the 
summations are being restricted of one unity in each further derivative. This 
will be clear in the last section, but you can conclude from the above calculations 
that: 
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' X <0 
a; < 1 aioPi(x); 



f'{x) ^{ x<2 anPnix - 1) + aoiPi(.T 



1) 
2) 



(32) 



a; < 3 ariPoix - 2) + ao2Pi{x 
x>3 



' X <Q 
X <1 

f"{x) = { x<2 





aooPo(a;); 
aoiPo(a; - 1) 
002^0 (a; - 2) 




(33) 



X < 3 
X > 3 

V 



Remark that Po can be replaced by x, but this is of no help at all. 

2.4 The forth convolution power 

From this point the recurrence is established: 

1. The convention for indexes follows the syntax of the language C++ for 
which the first positive (integer) index is zero, this way in a n x ti matrix 
the indexes run from ... n — 1. The matrix (a^ ); i, j = 0, • • • , n — 1 where 
n is the desired power; 

2. to use the entries in the last row of this matrix to define f{x) as combi- 
nations of translations of Pi{x), Pi{x — k), to the interval [k, k + 1]. 

3. The n — 1 derivatives of / are obtained, from the raw rt — 2 to the raw 
zero with increasing order of derivation and decreasing the number of 
summations of powers of P; in the same way as /. 

4. The starting point is the derivative of order n — 1 which comes from the 
line of order n — 1 of the alternate Pascal triangle. 

The forth power, /, of x ^Suy, now, be written, directly, without further 
integrations as an example before writing the algorithm formally. 

By Theorem [2l / has tree derivatives, two continuous derivatives, it is 3- 
splines. So the algorithm will start from the line of order 3 of the "alternate" 
Pascal triangle: 



aok 



Bin{3,k);k = 0...3; 



(34) 




aiiPo(l) + aoiPi(l); 
ai2Po(l) + ao2Pi(l); 



(35) 
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aiiPi(l) + aoiP2(l); 
ai2Pi(l) + 002^2(1); 




020 = aoo^o(l) 

021 = a20-P2(l) 

022 = a2i-Po(l) 

023 = 022^0(1) 

= aoofb(l) 
= 030^3(1); 

= a3iPo(l) + a2iPi(l) + aiiP2(l 
= 032^0(1) + a22Pi(l) + 012^2(1 



ao2P3(l); 



(36) 



(37) 



and 



X <0 

X < 1 

X <2 

X <3 

X < 4 

X > 4 





a30-P3(a:;); 
031^0(2; - 1) 
032^0(2; - 2) 
033^0(2; - 3) 




- 021^1(2; - 1) + aii^'2(-2; - 1) + 001^3(2; - 1) 

- 022^1(2; - 2) + ai2P2(-2; - 2) + 002^3(2; - 2) 
a23-Pi(2; ~ 3) + ai3P2(2; - 3) + ao3P3(2; - 3) 

(38) 

Observe that by Theorem [T] the derivatives of / can be traced back by 
difference of translations of the previous powers up to the last function-derivative 
as a linear combination of translations of x so we are able to right the equations 
of / and all its function-derivatives very easily. Better said, a program can do 
that for us, and the program can be download from a link at j4j. 



3 The n-th convolution power 

In the next theorem the indexes are separated with commas to avoid the 
confusion: "ri — 1, 1" for "n — 11". The same has to be done at the program. 



Theorem 4 (The n-convolution power) of the characteristic X[o,i] 



that 



The n-th convolution power of x ~ X[o.i]7 f = x"; is n ~ 1-b-splines such 



( 



X <Q 
X e [0, 1] 



71-1 

an-i,jPj{x - k);k ^ [x]] {1 < x < n) 
X > n 



where the coefficients a„_i.j are given by the equations: 



ai,o 



004 
ai,i 



ao,2 
ai,2 



fln-l.O an-1,2 



ao,n-i 

^n— l.n— 1 



0-o.k ~ Bin{n — 1, fc); fc ~ 0, . 



(39) 



l;Pm(n,p) = (-l)P("-i) (40) 
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ai.o = ao,o; (41) 

ai.i = ai,oA(l); (42) 

J >2;j<n- l;ai,j = ai,,_iFi(l); (43) 

a2,o = ao.ofo(l); (44) 

02,1 - 02,0^(1); (45) 

j >2;j< n;a2.j = 02 j_i + ai,j_iPi (1) + aoj-iP2(l); (46) 

03,0 = ao,o^o(l); (47) 

03.1-03,0^1(1); (48) 

j>2;j <n;a3,j= ^ bi,j-lPk{l) (49) 

i+fc=3 

(50) 

an-1,0 = ao,o-Po(l); (51) 
ari-1,1 = an-i,of'i(l); (52) 
J > 2; j < n - 1; a„_ij = ^ a,j_iPfc(l) (53) 

i+/c— n— 1 

The proof: 

By theorem\^ n-th convolution power, which is a \)-b-splines, has n—1 
function- derivatives, of which the first n — 2 are continuous, and the starting 
point is line of order n—1 of the alternate Pascal triangle defined by the numbers 
Bin{n — 1, fc). This is a consequence of two facts: 

1. In the combination of lemmas]^ and\^ we have a relation of the type of 
Stifcl rule, the defining relation of Pascal Triangle, modified to produce the 
alternate numbers Bin{n — 1, fc) = (~l)'^(fc~^)- So after n~l derivatives 
we are at line of order n ^ \ of the alternate Pascal triangle. 

2. Each new derivative is a linear combination of the previous convolution 
power by theorem [7] so the order of convolutions power will decrease to 
reach order one of derivation, the non- continuous function-derivative. 

This explains the first row of the (n — 1) x (n — l)-matrix b[ji]. The other rows 
are obtained by integrating the combination of translations of the Pj all along 
of the integers of the support [0,n] of f , and at each new integration the total 
variation of the previous integral is used as the initial value, with j being the 
order of the row, j ~ n — 1 — k where k is the order of the derivative. 



4 A computer program 

There is a program written in Calc, [3], which can be download from [4], 
and here you are the output of the program at the figure ([T]) for the choice of 
n = 30. This is the graphic of /, /', /" when / = x^'^. Time of this job: a few 
seconds. 

To have the graphics of /, /', /", of course, for n > 3, do the following: 
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30-th convolution power 



-0.5 
"dadosf" 
"dadosdl f" 
"dadosd2f" 
0.4 




Figure 1: graphics of convolution power 30, /,/'/" 



1. Start Calc, then you shall have a terminal of the calculator to make inputs. 
Calc is iterative, it is an interpreted language; 

2. Supposing you have the code at the current directory named to 

convolution_powers . calc 

execute at the terminal of Calc: 

read convolution_powers . calc 

"execute", means, "write the above at the terminal and hit enter". This 
will make Calc learn all the functions defined at convolution_powers . calc 
and you shall be able to execute any of them. There is one function called 
mainO which is written to do a big part of the job so you can have what 
is promised in the paper. 

3. Execute "main(n)" selecting the desired value of n at the terminal of Calc 
and you will see /, /', /". This will work out of the box, under Linux. 

You may experience some problems but then read the code and you shall 
see a suggestion to go around them. 
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4. If you want see the matrix, execute: 

• a = cria_matriz(n); either repeat the the value of n or put a new one. 

• imprimejnatriz (a,n) the parameter "a" has been created in the 
previous step, but you stiU need to put the same value of n as param- 
eter. This indicates that the program is not well planned... Perhaps 
somebody will produce something better from this draft. 

5. It is pointless to show all the n ~ 1 derivatives as it will be not possible 
to distinguish who's who for a large value of n. But, again, read the code 
and find a remark with the keyword "further derivatives" to see how to 
implement more derivatives. 



5 Final remarks 

As / = x" is a kernel, in the sense that it is positive and its integral is one, 
then it is integral preserving under convolution product. Well, this is exactly 
the reason why / has this property because x is a kernel. This permits to 
construct a partition of the unity from a set of characteristic functions associated 
with a partition of an interval with the desired class of differentiability. This 
can be done by making the convolution of / with each member of this set of 
characteristic functions and the methodology of sections 2,3,4 of paper can be 
of help to produce this rapidly. 

Given a kernel, in the above sense, we can transform it in a sequence which 
is an approximate identity, this is the terminology of [H chapter 6, definition 
6.31], the equation is = nf{nx). By the way, this is the expression to have 
a sharp kernel when constructing a partition of the unity, instead of using / 
directly. A simple translation will produce a kernel balanced around zero, and 
sometimes this is interesting to have. The program presented in this paper can 
be easily modified to make these applications. 

A simple consequence of theorem applied to the an n-splines, show 
that its derivative of order n ~ I will be a linear combination of characteristic 
functions of the elements of a partition of an interval of R (which can be the 
whole line). This permits a simple characterization of a quite general class of 
univariated splines. 



Theorem 5 (Characterization) spaces of splines 

Given a partition 11(7) of an interval lofR (which can be the whole line), 
and a kernel (f) being an n-splines of compact support 

1. let's consider 11(7) as the index set of the characteristic functions of the 
intervals in 11(7); 

2. the convolutions of (j) with the elements of 11(7) are linearly independent 
and generate a vector space of functions of class C"^^ whose dimension is 
card(n{I)), SplnAI); 
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3. the elements of Spln^cpil) is a space of n- splines. 

Observe that that there is a great generahty of objects which can be obtained 
this way, this diversity resulting from the possible shapes of </>. In case is the 
power of a characteristic function, the methodology of this paper applies to 
the construction of SplYi,4,{I)- You can have a glimpse of the power of this 
methodology looking at the graphic of [a bad example] 0] . 
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